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Fluorescence Intensity Multiple Distributions Analysis: 
Concurrent Determination of Diffusion Times and Molecular Brightness 

The present invention relates to a method for characterizing samples having fluores- 
cent pallidas, said method being based on the concurrent determination of diffusion times 
and specific brightnesses, or measures thereof. 

Fluorescence correlation spectroscopy (FCS) has proven to be a technique with single 
molecule sensitivity. Magde et al. (Phys, Rev. Lett. 29; 704 - 708, 1972) demonstrated the 
feasibility of detecting molecular number fluctuations by fluorescence correlation spectros- 
copy- Since then an increasing number of publications has appeared, aimed at improving 
the performance and accuracy of this technique. A progress was the implementation of con- 
focal detection optics (Koppel et al., Biophys. J. 16: 1315 - 1329, 1976; Rigler and Widen- 
pen, Bioscience 40: 180 - 183, 1990) and the use of silicon photon detectors (Rigler et al, 
Eur- Biophys. I 22: 169 - 175, 1993). This development pushed the detection limit below 
the single-molecule level (Rigler et ah, In; Fluorescence Spectroscopy: New Methods and 
Applications. O, S. Wolfbeis, editor. Springer.- Berlin. 13 - 24, 1993; Eigen and Rigler, 
Proc. Natl. Acad. Sci. USA §\ : 5740 - 5747, 1994; Brand et aL, J. Phys. Chem. 101 : 43 13 
-4321, 1997; EggelingetaJ., Proc. Natl Acad. Sci. USA 95; 1556- 1561, 1998). 

In recent publications comprising fluorescence fluctuation spectroscopy the attention 
has been diawn towards analyzing the histogram of the number of photon counts rather 
than the autocorrelation function (Qian and Elson, Biophys. I 57: 375 - 380, 1990; Fries et 
al, J. Phys. Chem. 102; 6601 - 6613, 1998; Chen ex al., Biophys. J. 77: 553 - 567, 1999; 
Kask et al, s Proc Natl. Acad. Sci. USA 96: 13756 - 13761, 1999).Lateiy, FCS has found a 
complement in the form of fluorescence intensity distribution analysis (FIDA). FIDA has 
been introduced as a method for analyzing mixtures of fluorescent particles. It is based on 
the detection of instantaneous photon emission rates from an open confocal volume. Cen- 
tral part of the method is the collection of photon count numbers, recorded in time intervals 
of fixed duration (time windows) and using this information to build up a count number 
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histogram. A theoretical probability distribution of photon count numbers is fitted against 
the obtained histogram yielding specific brightness values, q s as well as concentrations, c, 
for all different species in the sample. The term "specific brightness' 1 generally denotes the 
mean count rate of the deLeclor from light emitted by a particle of given species situated in 
a certain point in the sample, conventionally in the point where the value of the brightness 
profile function is unity. 

The historic predecessor of FIDA is FCS, which distinguishes different species on 
die basis of their characteristic diffusion times r, by analyzing the second order autocorre- 
lation function of light intensity, G(t) « (/(0) /(/)) - (/)* . Parameters which can be deter- 
mined by FCS (in addition to diffusion times x) are not, however, concentrations and spe- 
cific brightness values of different species separately, but products of the form cq\ While 
fluorescence intensity distribution analysis (FIDA) relies on a collection of instantaneous 
values of the fluctuating intensity, FCS analyzes the temporal characteristics of the fluc- 
tuations. Hence, the two methods represent complementary tools: FCS resolves compo- 
nents with different diffusion coefficients, while FIDA distinguishes the species according 
to their different values of specific molecular brightness. 

The present invention is aimed at developing a method that not only unifies the ad- 
vantages of both techniques, but rather yields a synergistic effect. 

According to the present invention, a new fluorescence fluctuation method is intro- 
duced that is based on the global analysis of a set of photon count number histograms, re- 
corded with preferably multiple widths of counting time intervals simultaneously. This 
Fluorescence Intensity Multiple Distributions Analysis (FIMDA) distinguishes fluorescent 
species on the basis of both, the specific brightness and the transitional diffusion time. 
The combined information, extracted from a single measurement, increases the read-out 
effectively by one dimension and thus breaks the individual limits of prior art techniques, 
i.e. FCS and FIDA. FIMDA can be widely applied for monitoring molecular interactions 
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including receptors and ligands or antibodies and antigens, which arc both of great rele- 
vance in the life sciences. 

The present invention describes the dependence of photon count number distributions 
on diffusion coefficients and applies it to a series of photon count number histograms cor- 
responding preferably to different widths of counting time intervals, in the following, 
sometimes the wording "experimentally determined distribution function" or "experimen- 
tally determined distribution" will be used instead of "histogram". As explained* the key is 
to analyze a set of distributions which is sensitive to the translational diffusion of particles. 
FCS detects the dynamics of particles^ because it compares the instantaneous intensities at 
time intervals separated by a certain delay. In order to make the distribution of photon 
count numbers sensitive to the temporal evolution of intensity one may alternatively - ac- 
cording to the present invention - choose to build a set of photon count number distribu- 
tions corresponding to different time windows. The choice of the time windows should 
preferably span a range comparable to the delay values used in FCS. 

In contrast to other two-dimensional FIDA techniques (Kask et al„ Biophys. J. 78, 
2000), which use two detectors, in FIMDA only a single detector is needed, though multi- 
ple detectors might be used. 

While (be ability of the method according to the present invention to determine spe- 
cific brightness values, diffusion times, and concentrations from mixtures is demonstrated 
on simulated data, its experimental utilization is shown by the determination of the binding 
constant of a protein-ligand interaction exemplifying its broad applicability in the life sci- 
ences. Since the values for all three parameters, i.e. brightness, diffusion and concentration, 
can be attributed to each species, the method according to the present invention is a more 
powerful tool of analysis than applying FCS and FIDA separately. The statistical accuracy 
of each extracted parameter determined according to the present invention is comparable to 
that in FCS or FIDA, often even better. 
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In the following, a modification of the theoiy of FID A is presented, which is a suit- 
able approximation for experimental purposes. In FIDA S a convenient representation of a 
pholon count number distribution P(n) is Us generating function, defined as 

a) 

The simple theory of FIDA assumes (i) that molecules are immobile during the 
counting time interval, and (ii) that the light flux from a molecule can be expressed as a 
product of a spatial brightness function B(t) (this is a function of spatial coordinates of the 
molecule characterizing the equipment) and a specific brightness q (characterizing a certain 
molecule species). Under these two assumptions, the distribution of the number of photon 
counts, emitted by molecules from a volume element dV, is double Poissonian and the cor- 
responding generating function reads 

W« s «P^*"* -fr)r -l)l> (2) 

where £ is the complex argument of the generating function, c is the concentration of mole- 
cules, and T is the width of the counting time interval. The representation is particularly 
convenient, because contributions from independent sources, like different volume ele- 
ments or species, are combined by simple multiplication of the contributing generating 
functions. The generating function of P(n) for a single species is 

*j. w G) = ^lcj{e {i ' m)r - l)dV], (3) 
while accounting for multiple species simply yields 



(4) 



The integral on ilic right hand side of Eq. 4 is calculated numerically, but instead of 
the three-dimensional integration over spatial coordinates, a one-dimensional integration 
coordinate x = ln|Z? 0 / B{r)} is introduced. The relationship between the brightness B and 
the coordinate x is therefore = J} 0 t< \ In FIDA it is .suitable to express the function 

dV Idx* which describes the brightness profile in onc-dimcnsional representation, by the 
formula: 



^ = 40 + + V> ( 5 ) 



Here ai, a 2 and a 3 arc empirical adjustment parameters granting for a sufficient flexi- 
bility la fit the measured histograms with high precision. However, a representation in the 
dV 

fcnii — = A 0 (x + a x x 2 + a 2 x 3 ) might also be used. The selection of coefficients Aq and Bo 

is nothing but the selection of the units of Vmd B. Usually, they are determined from the 
conditions 

J (6) 

J * CO 

So far, a simple version of the theory of FIDA has been described. For the purposes 
of FIMDA according to the present invention, one has to abandon the assumption that 
molecules are immobile during the counting interval. Surprisingly, according to the present 
invention one does not abandon Eq. 2, as well as the following equations, but qne will re- 
define the meaning of some variables instead, x is still a variable related to the spatial 
brightness profile, but now it characterizes the path of the molecule rather than its position. 
B is the spatial brightness averaged over the path rather than determined at a fixed position 
of the molecule. V is not the volume in space but dV/dx still expresses the probability that a 
molecule has a given value of x. If one would keep the original meaning of c and q s one 
would have to develop a theory predicting how Aq, a {i a 2 and <z 3 depend on the counting 
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' lime interval '/". However, according to the present invention another approach has been 
chosen. The normalization conditions have been kept (Eqs. 6 and 7) and it wis even found 
possible io apply a single selection of the values /f 0 , *vj and c/ 3 (or n .set of different time 
windows. The consequence of this selection is that in Rqs. 2-A c is an apparent concentra- 
tion (Capji) and is an apparent brightness (7^) which both depend on the width of the 
counting time interval T. 



In the following, a theory is presented predicting how and q m depend, on T- The 
case of single species is studied and the first and the second factorial cumulants of the dis- 
tribution corresponding to Eq. 3 are calculated. The factorial cumulanls are defined as 



(8) 



yielding: 

! ^=(«(«-d)-(«) j =c^v 2 

where normalization conditions given by Eqs. 6 and 7 have been used. (Note that Eqs. 9 
and 10 are in total agreement with Qian and Elson's formulae (Biophys. J. 57: 375 - 380, 
1990) derived under the assumptions i and ii.) From Eq. 9 one can conclude that 



(9) 
(10) 



c apF (?)g app (T) = (l), , * (11) 

where (/} s (n) r IT is the mean count rate, which does not depend on the choice of T. One 
shall proceed by employing the following relationship between the second cumulant of the 
count number distribution P(n\ T) and the autocorrelation function of fluorescence intensity 
G(0 = (/(0)/(0)-(/) 2 s 
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Introducing the nolaiioa 



t r 



one gets from Eqs. 12 and 10 

^ OX P 2 (7) - c app (Q) Up 2 (O)r(D. 
From Eqs. 1 1 and 14 one gets 

c m = -^^ 



(13) 



(14) 



(15) 
(16) 



As the concluding step in the theory underlying the present invention, one shall sub- 
stitute the expressions of G(t) fiom FCS into Eq. 13. If one ignores triplet trapping and 
studies pure diffusion, then c m (0) is the true concentration c , and q m (0) is the true spe- 
cific brightness q . Applying a Gaussian brightness function (Aragbn and Pecora, J. Chem. 
Phys. 64: 1 79 i - 1 803, 1 976), the autocorrelation function is 



Gw«) = cq 2 


Ml 
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(17) 
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denoting D as the diffusion coefficient and a, as ihc radial and a 7 , as the longitudinal dis- 
tance, where (he Gaussian profile has dropped e m times. The integrals in Eq. 13 yield the 
correction factor for translational diffusion 
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(18) 



where t = DT/a r 2 and /3 = or/ /a z 2 . For reasons explained below it is useful to calculate 
the first order terms in Eq. 18: 



1+ 
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+ 0(D 2 ). 



(19) 



However, from theoretical considerations as. well as from measurements it is known 
that simple physical models like Gaussian or else Gaussian-Lorentzian do not exactly rep- 
resent the actual brightness profile (Kask et al. 5 Proc. Natl. Acad. Sci. USA 96: 13756 - 
13761, 1999), Therefore, according to the present invention a modified Eq. 19 has been 
used and a fitting parameter a has been introduced, that preserves the first order tezms in 
Eq. 19: 



r«r(r). 
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(20) 



By matching the second order terms the Gaussian brightness profile would corre- 
spond to a = 2/3, but we rather choose a to be an empirical parameter,' which can be de- 
termined by the fitting procedure. From Eqs. 15 and 16 one can express the apparent pa- 
rameters of a pure diffusion process: 
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Solving Ihcsc equations with i aspect to f/^'and yields 



' . -, - (27) 

C ) CD = ^(7'). (28) 



Now, liaving solved the problems with diffusion and triplet transitions separately, 
one shall study the joint problem. Usually, the time scale of triplet transitions is much 
shorter than that of diffusion. Therefore, one is justified to replace c and q in Eqs. 21 and 
22 ^ ^ and $ jjn . This lets one combine Eqs. 21, 22, 27, and 28 to express c m and 

1w 515 



r c (T) = — - 

w r^re r (7T0+«)' (29 ) 



After having derived these expressions for c m and q m , the data simulations and the 
experiments should verify their validity. 

Taking into consideration the foregoing comments, the present invention relates in a 
general aspect to a method for characterizing samples having fluorescent particles, com- 
prising the steps of: 

- monitoring intensity fluctuations of fluorescence emitted by the particles in ac least one 
measurement volume by detecting sequences of photon counts by at least one photon 
detector, 



- determining from (he sequences of phoion counts intermediate statistical data com- 
prising at least iwo probability functions, f,( n i)'A( u ;i )»•••> of the number of photon 
counts, n M u^... a detected in different sets of counting lime intervals, 

- deuamining from said intermediate statistical data a distribution of particles as a func- 
tion of at least two arguments, wherein one argument is a (measure of) specific bright- 
ness of the particles and another argument is a (measure of) diffusion coefficient of the 
particles. 

In a prefeircd embodiment, the diffusion coefficient of a particle can be determined. 
It might however be also advantageous to use a measure which describes the diffusion of 
particles, e.g. diffusion time. The equivalent considerations apply to specific brightness, or 
a measure thereof. 

Generally speaking, the present invention therefore comprises the steps of 

(i) Optically sensing signals resulting from fluorescence emitted from fluorescent par- 
ticles like e.g. molecules present in the sample, and 

(ii) Electronically processing the optically sensed signals, wherein said processing 
comprises the step of concurrent determination of specific brightness and diffusion 
times of said particles based on probability functions of photon counts detected in 
different sets of counting time intervals. 

In a preferred embodiment, each set of counting time intervals consists of intervals of 
equal width while different probability functions P ?i (n,),^ (n 2 ) 3 ... correspond to counting 
time intervals of different widths r p T 2 In principle, at least two distributions with dif- 
ferent widths of the counting time interval axe nfceded in order to obtain information about 
the diffusion, but it might be preferred to use more than two distributions for analysis. 

In a further preferred embodiment, in each set of counting time intervals these inter- 



12 

" vals arc consecutive in lime. However, in each set of counting lime intervals these might, 
also overlap. 

It is particularly preferred that said distribution function of particles is determined by 
filling ihc experimentally determined probability functions /'(n,), A,(n 2 ).... by 
corresponding theoretical probability functions ^(fl,),/^!^)^ 

hi a preferred embodiment, said intermediate statistical data are processed applying 
inverse transforroation with regularizalion (ITR) or inverse transformation with constraints 
(ITC) or inverse transformation with both regularizalion and constraints (TTR.C). Because 
of statistical errors and limited size of measured data, inverse transformation as such is 
often an ill-posed mathematical problem, characterized by wild oscillations in its outcome. 
ITR, ITC and ITRC stabilize the mathematical problem by looking for a "regular" (e.g, a 
smooth) or constrained solution, for example by minimizing the sum of squared deviations 
of statistical data and a function of the solution itself, penalizing "irregular", usually ine- 
producible structures in the outcome, or values having no physical meaning. An example of 
constraining is disallowing negative values for concentrations (For the method of ITR, see 
e.g, W.H. Press et al Numerical recipes in C: the art of scientific computing, 2 nd edition, 
Cambridge University Press, 1992, p, 808, which is hereby incorporated by reference). 

In a preferred embodiment, the theoretical distributions /i(»i) a ii(n 2 ) a ... are 
calculated through their generating functions G Ha) (t)-Y^P(^)* In particular, said 

distribution function of particles can be determined by fitting the experimentally 
determined probability functions P 2 (« 2 ) 3 ... by corresponding theoretical probability 

functions iUM^O^X--, which are preferably calculated through* their generating 
functions G, w (fl = 

n 

In calculations of the theoretical distributions / > i (u | ),/ , 2 (n 2 )^. t the optical spatial 
brightness function B(t) might be accounted for by a separately determined relationship 
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between brightness B and volume dements dV, Depending on the position of a panicle in 
the measurement volume, the detected values of its specific brightness might vary, Espe- 
cially in cases of high polydispcrsity of the sample and/or low specific brighlncss differ- 
ences of the particles in the sample, il is preferred - in respect lo the accuracy of the meas- 
ured specific brightness values - lo take ihc spatial brightness function of the optical set-up 
into account. The relationship between the spatial brightness if and volume elements dV 
can be expressed through a variable x=\n{BJB) by a relationship 

= i4 0 (l + fl|* + fl2* ) x °*> where Bo is maximum brightness and Ao 7 a h a 2 and a$ are 
empirical parameters of the spatial brightness function. 

In a preferred embodiment, the generating function is calculated using the expression 
G^) = e5cp[j^^)|cf 3 r(5 (5 -^ ra ^ where c(q) is the density of particles with spe- 
cific brightness q, T is the length of the counting time interval, and JS(r) is the spatial 
brightness profile as a function of coordinates. 

In a further preferred embodiment, concentrations of particles are selected to be ap- 
proximately one or less moldcules per measurement volume. Experiments if performed at 
significantly lower concentrations than one particle per measurement volume would result 
in a slow acquisition of meaningful information because most of the data collection time is 
spent on waiting, i.e. with no particles in the measurement volume. In the other limit of an 
increasing high number of particles, photon count number distributions approach a Gaus- 
sian limit; therefore it is increasingly difficult to resolve contributions from different spe- 
cies. 

Preferably* either an avalanche photodiode or a photomultipiicr is Used as a photon- 
detector. It should be remarked that utilization of a single detector to collect distributions 
using a set of different widths of the counting time interval is favourite due to its fastness. 
However, in some instances it might be preferred to use at least two photon detectors, 
monitoring fluorescence of different wavelenghts or polarization. 
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The method according \o the present invention can be preferably applied in the con- 
ductance of homogeneous assays, i.e. mix-and-measurc assays with no need of separation 
slops. 

It can be e.g. applied in diagnostics, high throughput drug screening, biomcdicinc, 
optimization of properties of molecules and identification of specific cell populations. 

According to the present invention, confocal techniques arc particularly suited to 
monitor fluctuating intensity of fluorescence. They may be applied - as outlined above - to 
a wide field of applications, such as biomedicine, etc, The conjugate focal (confocal) tech- 
nique is based on using a point source of light sharply focused to a diffraction-limited spot 
on the sample. The emitted light is viewed through a spatial filter (pinhole) that isolates the 
viewing area to that exactly coincident with the illuminating spot. Thus, the illumination 
and detection apertures are optically conjugated with each other. Light originating from 
focal planes other than that of the objective lens is rejected, which effectively provides a 
very small depth of field. Therefore, in a particular preferred embodiment of the present 
invention, a confocal microscope is used for monitoring the intensity of fluorescence. In 
order to achieve a high signal-to-noise-ratio, it is useful to monitor fluorescence using an 
apparatus that comprises; a radiation source (12) for providing excitation radiation (14), an 
objective (22) for focussing the excitation radiation (14) into a measurement volume (26), a 
detector (42) for detecting emission radiation (30) that stems from the measurement vol- 
ume (26), and an opaque means (44) positioned in the pathway (32) of the emission radiar- 
tion (30) or excitation radiation (14) for erasing the central part of the emission Radiation 
(30) or excitation radiation (14). ft might be particularly preferred to use an optical set-up 
described in detail in Figure 7. 

The advantage of FMDA (present invention) over FCS (prior ait) is that the method 
yields genuine concentrations of components in the sample, instead of the products of con- 
centration and brightness squared in FCS. Only the independent determination of at least 
one of the two specific brightness values enables FCS to determine two concentrations un- 



15 



ambiguously. However, inexperienced users of FCS often silently assume equal specific 
brightness when resolving two components. This assumption can cause significantly biased 
results, FIMDA brings this issue to the focus of analysis. 



Another advantage of the presented method is its versatility. If FCS (prior art) or 
FIDA (prior art) fail to detect a particular readout upon a biochemical reaction FIMDA 
(present invention) might be able to succeed. The biochemical reaction is not necessarily 
limited to the binding of two components, but can be any chemical reaction of interest The 
possibility of using only one detector for recording two physical characteristics in a single 
measurement makes FIMDA a very efficient method of analysis which saves precious as- 
say development time. 

Because two read-out parameters are recorded in one measurement preferably with 
one detector, the present invention is an efficient and reliable method of analysis. The in- 
troduction of the generation function facilitates data evaluation and makes the method a 
fast analytical tool. FIMDA is more universal applicable to biochemical assays than FCS 
and FIDA applied separately, because each assay component can be distinguished by both 
diffusion time and brightness. The statistical accuracy of determining the diffusion time in 
FIMDA is better or at least as good as in FCS. 

The invention is illustrated by the following figures and examples, which are not in- 
tended to limit the scope of the invention. 



Figure 1 illustrates the count number distributions and fits of a 3.8 nM Cy5 solution 
recorded simultaneously at different time windows T. The weighted residuals for the dif- 
ferent time windows are shown in the lower part of the figure. 



Figure 2 illustrates the fitting results of simulated data for a mixture of 3 compo- 
nents. The simulated brightness (in kHz) and diffusion time (in (is) values for the compo- 
nents are: (30 kHz, 192 us); (120 kHz, 192 us); (120 kHz, 64 |is). The contributions to the 
total intensity are 10.8 kHz f 20.4 kHz, and 14,4 kHz, respectively. The graph presents the 
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results of FIMDA from 20 independent realizations of annul atiuns, each correspond uig to 
an experiment of 60 s duration, 

Figure 3 illustrates the binding of pTyr-Val-AsivV^l-Ly5(Cy5) to SH2. The solid 
curve results from a hyperbolic lit, yielding a binding constant of A' 0 = 1.54 ± O.J 4 jxM. 

Figure 4 illustrates the experimental set-up used in Experiment 2. Radiation emitted 
by a laser passes an OD filter and reaches a dichroie mirror which reflects the radiation 
towards an objective having its focus within the sample under study. Fluorescence emitted 
from the sample passes the objective and reaches the dichroie mirror which is transparent 
for the fluorescent emission. After passing a bandpass filter and a pinhole, the emission 
reaches an avalanche photo-diode used as pan of the detector. By means of a photon 
counting unit and a computer, specific brightness and diffusion can be determined accord- 
ing to the present invention. 

Figure S shows 10 count number distributions with time windows 40, 60, 120, 200, 
400, 600 4 800, 1200, 1600 and 2000 \ls from a 0.8 nM Cy 5 solution. 

Figure 6 shows the calculated apparent specific brightness of the dye as a function of 
counting time interval; evaluated by FID A. 

Reference is now made to Figure 7 which shows one embodiment of an apparatus 
adapted for use in performing the method according to the present invention. Apparatus 10 
comprises a laser 12 which serves as a light source for illuminating the sample by a bundle 
of coherent monochromatic excitation radiation 14, Excitation radiation 14 is paralleled by 
a lens 16 and reaches a dichroie mirror 20. Preferably, the angle between, the optical axes 
18 and the dichroie mirror 20 is 45°. The dichroie mirror 20 reflects the excitation radiation 
14 in direction of an objective lens 22 having its focus 24 within a sample volume 26. 
Sample volume 26 and objective lens 22 are preferably separated from each other by a 
transparent cover glass 28, e,g. by the bottom of a commercially available micro-titer plate 
which houses the sample. The sample preferably includes fluorescently labeled molecules 
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or oilier particles. Due to excitation by an appropriate excitation radiation 14 ? ihc molecules 
or other particles present in the sample emit radiation 30. Emission radiation 30 passes the 
objective lens 22 and reaches the dichroic mirror 20 which is transparent for emission ra- 
diation 30. Thereafter, emission radiation passes a filter 34 and a collimator lens 36 on the 
optical axes 32. A pinhole 38 is situated in ihc focus of collimator lens 36. Emission radia- 
tion 30 passing the pinhole 3S reaches a further lens 40 and, thereafter, is detected by the 
photo-detector 42. Within the pathway of emission radiation 30, in particular between di- 
chroic mirror 20 and photo-detector 42, an opaque means 44 is provided through which a 
central part of the emission radiation 30 cannot pass. This central part of the emission ra- 
diation 30 stems from areas on the optical axes 32 in front of or behind the focus 24 of the 
excitation radiation 14. Only emission radiation 30 that stems from the focus 24 or its di- 
rect neighborhood passes the pinhole 38 and reaches photo-detector 42. Instead of placing 
an opaque means 44 within the pathway of emission radiation 30, the pathway of excitation 
radiation 14 is also suitable for positioning an opaque means 44. In particular, an opaque 
means 44 can be positioned between laser 12 and dichroic mirror 20. The use of an opaque 
means 44 in the rixethod according to the present invention as described in detail herein 
improves the signal-to-noise rati o. 
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EXPERIMENT 1 
MATERIALS AND METHODS 
Experimental set-up 

The central optical pari of the present F1MDA expcrimcnl is a confocal microscope 
as it is used in fluorescence correlation spectroscopy (KLoppel cL q]., Biophys, J. 16": 1315 - 
1329, 1976: Rigler ct aL, Eur. Biophys. J. 22: 169 - 175, 1993), For the excitation of fluo- 
rescence, a beam from a continuous wave laser is attenuated to -800 jiW by neutral densixy 
filters, passes a beam expander, and is directed to the microscope objective (UApo/340* 
40x, N.A. 1.15, Olympus Optical Co. Ltd, Tokyo, Japan) by a dichroic mirror. Fluores- 
cence is collected by the same objective through the dichroic mirror, a spectral band-pass 
filter, and is focused to a confocal pinhole, which selves to reject the out-of-focus light 
The light, which passes the pinhole, is detected by a silicon photon counting module 
(SPCM-AQ-131, EG&G Optoelectronics, Vaudreuil, Canada). An electronic counter, con- 
structed at EVOTEC as a computer plug-in card, collects the TTL pulses from the detector 
continuously and calculates the count number histograms for all preselected widths of time 
windows (40, 60, 120, 200, 400, 600 3 800, 1200, 1600, 2000 ps) in real-time from the 
32 MB onboard buffer. By feeding the detector outputs to a correlator, FCS measurements 
can be performed in parallel with FIMDA experiments. 

In order to satisfy the spectral needs of the various fluorophores used in this study, 
different lasers and spectral band-pass filters were employed. For Cy5 (Aroersham Pharma- 
cia Biotech, Bucks, UK) conjugated bio-molecules an arrangement of a red laser diode 
(Crystal GmbH, Berlin, Germany; 635 nm) and a band-pass filter with a central wavelength 
of 670 nm (670DF40, Omega Optical, Biattleboto* VI) was used. In case of TAMRA (5- 
carboxytetramethyirhodamme) labeled molecules this was an arrangement of a frequency 
doubled Nd-YAG laser (jiGrccn 4601; Uniphase, San Jose, CA; 532 tun) and a 59GDF6Q 
filter. 

The effective dimensions of the illuminated volume were calibrated indirectly, using 
FCS on small dye molecules (TAMRA, Cy5) with known diffusion coefficients. The auto- 
correlation functions of diffusion were fitted to Eq. 17, i.e. assuming a 3-dtmensionaI 
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Gaussian intensity profile. The exact determination of the dimensions and profile would be 
very complex, since they are affected by both, the size of the laser beam as well as the size 
of the confocal pinhole, However, in most cases the knowledge of the exact dimensions is 
not necessary. 

The focal beam radius was adjusted to -0.75 jwn by selecting an appropriate expan- 
sion factor of the original laser beam, resulting in a mean translational diffusion time of 
360 us for the free dye Cy5, This diffusion can be clearly observed when rising the time 
windows from 40 ps to 2 ms. As can be seen in Fig. 1, the selected count number distribu- 
tions of a 3.8 xM CyS solution differ considerably. However, die major differences be- 
tween the distributions are due to the varying mean count number in different time win- 
dows used. Diffusion of fluorescent molecules causes only small but significant modifica- 
tions to the shape of each distribution. 

The levels of background count rate are determined in a separate experiment on 
bidistilicd water and amount usually to 0.5 kHz. The main contributor to this non- 
fluctuating background' light intensity is Raman scattering from water. 

Data simulations 

Real samples, comprising a mixture of molecules, which express deliberately chosen 
parameters (brightness values and diffusion coefficients), are difficult to prepare. There- 
fore, some evaluations of the new method were performed using simulated data. A number 
of sets of histograms for FIMDA, HDA, and correlation functions for FCS have been 
simulated according to the following algorithm. In a closed rectangular reservoir, a given 
number of molecules is initially randomly distributed over a high number (typically 
360*360*720) of discrete spatial grid points. Each molecule is subject to' consequent diffu- 
sion simulation and jumps randomly by one grid unit either in x-, y- 4 or z-direction with a 
frequency corresponding to a given diffusion coefficient The "focus" is located in the 
center of the reservoir, and the brightness distribution is assumed to be Gaussian in all 
three dimensions. When calculating the brightness integral from a molecule over a given 
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set of time intervals, Uic molecule can be randomly trapped* as well as released from the 
triplet excited slate (where it is dark). Now one can calculate an array of brightness inte- 
grals over basic time intervals of a given width (e.g. 5 jxs) describing the evolution of the 
mixture. The brightness integrals arc then converted into photon count numbers by gener- 
ating n random l'oisson number with the corresponding average. This step also accounts tor 
the noise introduced by the detector because the random number generator is used not only 
for driving random motion of molecules but also for simulating random numbers of de- 
tected photons at given light intensities. The random count numbers obtained were subse- 
quently used to calculate histograms for FIMDA and FIDA as well as the correlation func- 
tion for FCS. 

Due to the finite size of the simulation reservoir, some distortions of the correlation 
function (Le. s deviations from Eq. 17) can be expected. The distortions axe in fact below 
the statistical noise level. Therefore, one can consider the simulations to be an adequate 
tool for estimating statistical errors of the extracted parameters. For this purpose, typically 
30 realizations of experiments with a given set of molecular parameters were simulated, 
from which the standard deviations and the coefficients of variation (CV) as the ratio of 

standard deviation to mean value were calculated. 

\ ' 

Fitting 

A series of simultaneously measured or simulated distributions is globally fitted us- 
ing a Marquardt algorithm. The fitting program is a modest modification of the program 
designed for FIDA (Kask et al. f Proc. Natl Acad Scl USA 96: 13756 - 13761, 1999). 
Theoretical distributions are calculated using exactly the same algorithm as in FIDA, ex- 
cept that each species has an individual apparent concentration and an apparent brightness 
at each time window, calculated according to Eqs. 29, All parameters not assigned to spe- 
cies but rather to the equipment (Le.. A Q> a u a 2l .and a* from Eq. 5 and a from Eq. 20) are 
usually determined beforehand from separate adjustment experiments on pure dye solu- 
tions. 
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Biochemical System 

This Grh2 (SII2)-phoxphopq)ikk interaction 

Recem niiiikimor research Iras been focused on tyrosine kinase growth factor rcccj> 
tors (Levhzki, Eur, J. Biochcm. 226: 1- 13, 1994; Alcssandro ei al„ Curr. Top. Microbiol. 
Immunol. 213: 167 - 188, 1996; Furct ct al... J. Med Chem. 41: 3442 - 3449 : 1998). A 
critical link in the signal transduction pathway of (his receptor is the interaction of its phos- 
photyrosine residue (pTyr) with the Src-homology 2 (SH2) domain of the adapter protein 
Grb2 (growth factor receptor-binding protein 2). For the recognition, a minimal peptide 
sequence of the receptor (pTyr-Val-Asn) is sufficient (MiiUer et aL, J, BioL Chcm. 271: 
16500 - 16505, 1996; Gram ct aL, Eur. J. Biochern. 246: 633 - 637, 1997; Furct et al., J. 
Med. Chem. 41: 3442 - 3449 7 1998). The binding partner of this peptide motive, the SH2 
domain of Grb2, can fold into a functional protein module independent of neighboring se- 
quences (Booker et aL, Nature, 358; 684 - 687, 1 992; Overduin et aL, Cell 70: 679 - 704, 
1992). Therefore, as a model system, the bare SH2 domain (14.3 kDa) has been chosen to 
interact with a fluorescently labeled phosphopeptide (pTyr-VaKAsn-VaKLys(CyS)) (1387 
Da). 

The SH2 domain of Qib2 was prepared as described elsewhere (Lowenstein et al., 
Cell 70:431 - 442, 1992; Baumann et aL, Eur. I ImmunoL 24: 1799 - 1807. 1994; Mtiller 
et aL, J. BioL Chem. 271: 16500 - 16505, 1996). The phosphopeptide was synthesized 
using manual Fmoc solid phase chemistry and labeled with Cy5-NHS via a Lysine residue. 
An additional Valine was introduced to minimize possible interactions of the dye with the 
main recognition motive pTyr. The final compound, pTyr-Val-Asn-Val*Lys(Cy5) was 
characterized by mass spectrometry (LC/MS, and MALDI/TOF), UV/VIS, and fluores- 
cence spectroscopy. 

RESULTS 

Data simulations and test experiments 

At first, a series of measurements on a 1 nM TAMRA solution was performed col- 
lecting data in parallel, for FIMDA as well as for FCS. This scries of experiments, with 
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Juration of 2 s each* was repeated in simulations using .similar molecular parameters, Tl^o 
purpose of these experiments was to verify whether simulations arc <i reasonable model of 
real experiments., in particular whether data simulations arc u reasonable means of predict- 
ing statistical errors of eslimatcd parameters, The coefficients of variation of the parame- 
ters extracted from simulated data indeed coincide wiih the results of the real experiment, 
as can be seen in Tab. 1. 

Another scries of test experiments was repeated in a significantly shorter time do- 
main with the goal to compare FIMDA and FCS in their ability to estimate parameters of 
the triplet component A set of counting time intervals of 2, 4 3 8, 16, 32, 64, 128, 256, 512, 
and 1024 us was selected for this purpose. The duration of these experiments was 1 6 s. The 
results, presented in Tab. 2 4 indicate that the values for the triplet parameters estimated by 
FIMDA have similar dependence on the excitation intensity to the FCS results. The 
FIMDA results are slightly biased and have higher CV values compared to FCS, since the 
estimation of triplet parameters in FIMDA is indirect, because the shortest time window (2 
us) is equal to the triplet lifetime. However, the main purpose of the triplet correction in the 
model is not to determine the triplet parameters, but to improve the quality of the fit and to 
remove a source of bias in the, brightness and diffusion parameters. 

Out of curiosity histograms for FIMDA were simulated for three-component analy- 
sis. Two of the components had equal brightness values (120 kHz) ? and another pair had 
equal diffusion times (192 us). Due to the larger number of free parameters, the simulated 
duration of experiments was increased to 60 s, so that the variations of fitted parameters 
stayed in reasonable limits. In this test, all parameters were subject to fitting- The results 
are presented in Fig. 2 as vertical bars in a plane with brightness and diffusion time as x-y 
coordinates, and the ordinate displaying the contribution to the intensity, i.e r , the product of 
concentration and brightness. The three components are clearly resolved, since the scatter 
in the location of individual bars is much smaller than the distance between the groups, 
which correspond to different components. 

Note that with FIDA alone, the components with equal brightness cannot be resolved, 
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while with FCS alone, the components with equal diffusion lime remain unresolved. 
Biochemical .system 

Tiic experimental utilization of the new method will be demonstrated by ibc determi- 
nation of the binding constant of the above introduced Grb2 (SI12;-phosphopcplide inter- 
action. For ihis purpose <\ titration experiment was carried out keeping the pTyr-Val-Asn- 
Val-Lys(Cy5) concentration constant at 0.4 nM, while SH2 was subject to titration (0.01, 
0.03, 0,1, 0.3, 1, 3, 10, 30, 100, and 130 jiM). All experiments were performed under iden- 
tical conditions, i.e., the same buffer (sterile filtered water, 50 mM Na-phosphatc buffer pH 
7-8, 50 mM NaCl, and 0.05% Pluronic; T=20°C), and a data acquisition time of 30 s per 
measurement, repeated 30 times per sample* In each single measurement the same set of 1 0 
different time windows was used (40, 60, 120, 200, 400, 600, 800, 1200, 1600, 2000 \is) 
resulting in 10 different photon count number histograms, which were globally fitted. 

As the first stop, the diffusion time Ti = 407±6jis and the specific brightness 
q\ - 3 1.7 ± 0.3 kHz were determined from a single component analysis applied to the pure 
conjugate solution. Addition of excess SH2 (130 jiM) to 0.4 nM conjugate resulted in a 
sample with the majority of #ie conjugate bound to SH2. The complex was characterized 
both by a longer diffusion time and a higher specific brightness compared to the free con- 
jugate. This mixture was then analyzed by all three methods (FIMDA, F1DA, and FCS) 
using a two-component fit with i { and/or q { fixed, depending on the method The results of 
this step of analysis are presented in Tab- 3. It can be seen that all methods yield similar 
values of parameters for the complex. The corresponding CV values were again determined 
by two independent methods, i.e., from the statistical analysis of the results of a series of 30 
measurements and from simulations. The two estimates of the statistical errors agree rea- 
sonably well and the CV values corresponding to different methods are similar, with the 
exception of FID A, which has difficulties due to the small (30%) difference in specific 
brightness of the two components. 

As the next step of the present studies, a sample with 3 p.M SH2 was analyzed. This 
particular concentration was chosen to achieve a mixture of approximately equal proper- 
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tions of complex and free conjugate. Since it is rather difficult to resolve components with 
only a twofold difference in diffusion coefficient and even smaller difference in specific 
brightness, here also the diffusion lime and brightness of the complex were fixed to the 
values of Tab. 3. With the molecular parameters fixed, ihc concentrations were reliably 
determined by all methods. The results of this sicp of analysis arc summarized in Tab. 4. 

i 

In the same manner, the whole series of SH2 concentrations was fitted. Fig. 3 shows 
the calculated traction bound (c^icx / (complex + ^coni^ic)) for FIMDA with the solid curve 
resulting from a hyperbolic fit, that yielded a binding constant for the SH2~phosphopeplide 
interaction of K 0 = 1.54 ± 0.14 jjlM> Comparable binding curves were also obtained by FCS 
and FIDA (data not shown), with K D values of 2*16 ±0.19 (-iM and 1,60 ± 0.19 uM 3 re- 
spectively. 

The data of Fig. 3 demonstrate that FIMDA is a suitable method for monitoring the 
formation of a molecular complex. FCS and F1DA experiments yielded similar Ko values 
for this particular SH2-phosphopeptide interaction. In the literature the affinity is reported 
to vary by several orders of magnitude, depending on the peptide sequence (Muller et aL, J. 
BioL Chora. 271: 16500 - 16505, 1996; Gram et aL, Eur. J. Biochern. 246: 633 - 637, 
1997; Furet et aL, J_ Med, Chem. 41 - 3442 - 3449, 1998). High affinities are in the range of 
Ko 888 10 - 100 nM- However, with a lysine (and Cy5 attached to it) at the +4 position of the 
phosphopeptide (defining p-Thr as the 0 position with c +» continuing on the C and on 
the N terminus) the affinity decreases to the micromolar range. This result agrees well with 
the importance of lipophilic groups attached to 'appropriate' positions on the C-terminus, 
increasing the binding constant to the SH2-domain (Furet et aL, J. Med. Chem. 41 : 3442 - 
3449, 1998), For example Val (at position pTyr+3) is making van der Waals contact with a 
large hydrophobic area on the SH2-domain. 

One of the surprising results of this study according to the present invention is that in 
each of the experiments, the statistical accuracy of the diffusion time estimated by FIMDA 
is as good as or even better than that estimated by FCS. This is a counter-incuitive result 
because FCS is directly focused on fitting a diffusion-dependent correlation function G(t), 
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while in FJMDA the diffusion lime is estimated only indirectly, namely tluough the de- 
pendence uflhc apparent brightness on the width of (he lime window. 

A further observation in this respect is that the CV values of Ihc diffusion times arc 
in general higher than those for the brightness values. This also holds true for the theoreti- 
cal simulations and therefore reflects an effect rooting in the measuring principle. The phe- 
nomenon can be explained qualitatively by the different ways how these quantities are de- 
termined. For simplicity, one may imagine an observation volume with a constant bright- 
ness profile B(r) inside. In this case, one only needs to measure the average count rate of a 
molecule that enters the volume to determine its specific brightness. This requires the de- 
tection of many photons per given time interval but can in principle be achieved from a 
single passage. On the other hand, for estimating the diffusion time, one has to determine 
the mean duration of the diffusion driven passage, which inevitably requires averaging over 
many events, even though many photons may be detected each time. Therefore, in an ex- 
periment of fixed duration, the specific brightness of a molecule can in principle be deter- 
mined with a higher accuracy than its diffusion time. 

EXPERIMENT 2 

Figure 4 shows'the experimental set-up used in this study. As an example for a typi- 
cal measurement Figure 5 shows 10 count number distributions with time windows 40, 60, 
120 3 200, 400, 600, 800, 1200, 1600, 2000 (is from a 0.8 nM Cy 5 solution. Figure 6 
shows the calculated apparent specific brightness of the dye as a function of counting time 
interval (time window), evaluated by FID A. Since the diffusion time of the dye is long 
(compared to 40 jia) it spends only a short time in the confocal volume and its apparent 
brightness decreases. From this information the diffusion time can be calculated from a 
global fit. 



The viability of the present invention was tested by measurements characterizing a ii- 
gand-prolein interaction. The minimal peptide sequence pTyr-Val-Asn (originating from 
tyrosine kinase growth factor receptors) binds to the SH2-domain of the adapter protein 
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Grb2. The conjugate pTyr-VaI-Asn-Val-tys(Cy5) was synthesized as binding partner for 
SH2. All measurements were performed under identical conditions (e.g. same buffer: ster- 
ile filler water, 50 mM Na-phosphatc buffer pN 7.8, 50 mM NaCJ and 0,05 % Pluronic; T 
" 20 °C). To obtain the Ko value for the ligand protein interaction a titration experiment 
was carried out. The pTyr-Val-Asn-Val-Lys(Cy5) concentration was kept constant at 0.8 
nM, while SH2 was subject to titration (0.01, 0.03, 0.1, 0.3, i, 3> 10, 30 r 100, and 130 ^M). 
Each sample was measured 30 times keeping the acquisition time at 30 seconds per meas- 
urement FCS (prior art), FIDA (prior art) and FIMDA (present invention) data were col- 
lected simultaneously in order to compare all three methods. Table 5 summarizes the re- 
sults. It becomes clear that FIMDA is able to resolve both, the specific brightness and the 
diffusion time of both species (ligand and complex). Therefore, FIMDA is a more versatile 
tool than FCS or FIDA and will save assay development time. Comparing the standard 
deviations of all methods in determining brightness and diffusion time FIMDA is better or 
at bast as good as FCS or FIDA. 
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TABLES 

Tabic 1 Comparison of coefficients of variation of estimated parameters from .scries 
of experimental and simulated histograms by FIMDA, and correlation functions by FCS. 

Parameter Mean value CV (%) 



HMD A .■ FCS 







Experimental 
data 


Simulated 
data 


Experimental 
data 


Simulated 
data 


Brightness q (kH2) 


115 


2.5 


2.3 


2.6 


2.5 


Concentration c 


0.73 


3.1 


3.8 


3.8 


4.0 


(molecules per con- 












focal volume) 












Diffusion time-r(ns) 


287 


5.4 


4.3 


7.3 


5.2 



Table 2 Triplet parameters, estimated from a series of experiments on 1 nM TAMRA 
solution by FCS and FIMDA at two different excitation intensities. Excitation wavelength 
532 nm, duration 16 s, time windows 2, 4, 8, 16, 32, 64, 128, 256, 512, and 1024 us. 



Peak ex- 






FCS 






FIMDA 




citation 


Triplet 


CV(«/ 


6) Triplet 


CV(%) 


Triplet 


CV(%) 


Triplet 


CV (%) 


intensity, 


lifetime 




popula- 




lifetime 




popula- 




(kW/cm 2 ) 


0*0 




tion 




(Ms). 




tion 




118 


1.98 


3.9 


0.182 


2.7 


3.12 


7.6 


0.137 


3.1 


187 


1.75 


3.7 


0.235 


1.6 


2.59 


3.8 


0,183 


2.0 
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Tabic 3 Comparison of estimated parameters and iheir coefficients of variation at a. 
high receptor concentration (130 p,M).A scries of 30 experiments of 30 s durution each 
was evaluated by FIMDA* FfDA and KJS. Brightness (in FIMDA and PI I) A) and diffusion 
lime (in FIMDA and FCS) of the free conjugate were independently determined and fixed 
to 3 L7 kHz and 407 jis, respectively, ui this analysis. 



Parameter 


Method 


Mean value from 
experiment 


CV (%) 


CV (%) 
from simulations 


Cconjuguc 


F1MDA 


0.132' 


43* 


82* 


(molecules per 


F1DA 


0.196* 


76* 


71* 


confocal volume) 


FCS 


0.052* 


99* 


120* 


^complex 


FIMDA 


0.618 


9.6 


8.0 


(molecules per . 


FJDA 


0.555 


26.7 


14.9 


confocal volume) 


FCS; 


0.710 *" 


7.9 


12.6 




FIMDA 


39.5 


2.2 


2.3 


loomplex (kHz) 


FIDA 


38.4 


5.5 


3.9 




FCS 


36.4 


3.7 


3.4 




FIMDA 


0.913 


6.9 


4.6 


fcomplex(ms) 


FCS 


0.898 


5.4 


7.2 



out of 30 realizations, in 5 to 10 cases zero conjugate concentration was yielded by the 
fitting program (negative values are disallowed). This indicates that in this particular ex- 
ample the conjugate concentration could not be properly determined. However, for further 
data analysis only q mmr k X and T^mpic* are needed. 
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Tabic 4. Comparison of the estimated concentrations at an intermediate receptor 
concentration (3 pM). In addition to the brightness and the diffusion lime of the free conju- 
gate, also Lhc brightness and/or the diffusion time of the complex were fixed here to values 
shown in Tab. 3 



Parameter 


Method 


Mean value from 
experiment 


CV (%) 


CV (%) 
from simulations 


c conju&jie 


F1MDA 


0.328 


11.7 


4.9 


(molecules per 


FIDA 


0.311 


16.1 


4.9 


confocal volume) 


FCS 


0.303 


8.4 


8.0 


c complex 


FIMDA 


0.437 


8.2 


3.4 


(molecules per 
confocal volume) 


FIDA 


0.455 


10.0 


3.7 




FCS 


0.467 


7.8 


8.8 
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Table 5: Comparison of FIMHA (present invenrion) willi prior art FCS and l-HM 
with respect in the determination oi'q. t mid Kg 



[kHz] [kHz] 

F1MDA 31.7±0.32 39.6±0.9 
FCS 

FIDA 32.1*0.31 43.3±2.3 



[ms] [ms] " friM) 

0.402±0.006 1.01*0.08 L68± 0.27 
0.402±0.013 0.867±0.057 2.26±0.28 

1.70±0.29 



